THE LOCAL TRACE FUNCTION FOR SUPER- WAVELETS 



DORIN ERVIN DUTKAY 



Abstract. We define an affine structure on {R)(B...(BL^ (R) and, following 
some ideas developed in IDutl| . we construct a local trace function for this 
situation. This trace function is a complete invariant for a shift invariant 
subspace and it has a variety of properties which maJte it easily computable. 
The local trace is then used to give a characterization of super-wavelets and to 
analyze their multiplicity function, dimension function and spectral function. 
The "nx" oversampling result of Chui and Shi CS is refined to produce 
super-wavelets. 



Contents 



1. 


IntroductioiJ 


2. 


Pcrniutativc wavelet reoresentationsl 


3. 


Shift invariant subsoaced 


4. 


The local trace function] 


5. 


A characterization of suDcr-wavelctd 


6. 


The SDCctral function and the dimension functioni 


Referenccsl 



1 

4 
5 
8 

12 
18 
21 



1. Introduction 

The wavelet theory involves the study of an affine structure existent on a Hilbert 
space H, given by two unitary operators, T the translation, and U the dilation that 
satisfy the commutation relation: 

where iV > 2 is an integer called the scale. For the classical wavelet theory, the 
Hilbert space is (K), the translation is 

ro/(x) = /(x-i), (xeRJeLh 

and the dilation operator is 



where iV > 2 is an integer. 

Consider an affine structure U, T on a Hilbert space H. A wavelet is a set 
= {tpi, V'l} such that the affine system 
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is an orthonormal basis for H. Often the requirement is weaker, and one looks for 
the afSne system to be only a frame. We recall that a family {ei)iei of vectors in 
a Hilbert space H is called a frame with constants A > and S > if 

A\\lf<Y,\{f\e.)\^<B\\f\\\, ifeH). 
iei 

When only the second inequality is satisfied, the family is called Bessel. When 
A = B — 1, the frame is called a normalized tight frame. 

Each orthonormal wavelet has a generalized multiresolution attached to it. 

A generalized multiresolution analysis (GMRA) is a collection {Vj)j^z of sub- 
spaces of H that satisfy: 

(i) Vj C Vj+i for aU j. 

(ii) UVj = V,_i for aU j. 

(iii) UVj is dense in H and nVj = {0}. 

(iv) Vq is invariant under Tk for all fc S Z. 

If, in addition, there is a vector £ Vo such that {T'^ip | fc e Z} is an orthonormal 
basis for Vq, then {Vj)j is called a multiresolution analysis (MRA) and ip is called a 
scaling vector for this MRA. For more information on GMRA's we refer to |BMM| . 
The GMRA associated to a wavelet ^E* is defined by 

Vj {U"'T''tl; \m>-j,ke Z}, (.? G Z). 

If the GMRA associated to is actually a MRA, then is called a MRA wavelet. 

As it is shown in |HLj and |Dut2j . it is possible to construct an affine structure 
and wavelets on the larger Hilbert space (M) © ... (R), i.e. "super- wavelets" . 
One way to do this is by considering the operators: 

T:=ro® ...©To, C/:=C/o©...©C/o, 

and as shown in |HL| . wavelets can be constructed for this affine structure, but 
none of them are MRA wavelets. 

The affine structure described in |Dut2| takes into consideration the cycles of the 
map z 1-^ . A set {zi, Zp] of points on the unit circle T := {z G C | |z| = 1} is 
called a cycle if the points are distinct and z^ — Z2, Zp_i = Zp, z^ = zi. Then 
the translation is defined by 

T(/i, /p) = (ziTo/i, ...,ZpTofp), (/i, ...Jp e (M)), 

and the dilation is defined by 

uih, fp) = ({/0/2, C/0/3, Uofp, Uoh), (/i, fp e (R)). 

Note that the dilation operator permutes cyclically the components /i, fp. Then 
one can do finite direct sums of these representations for different cycles to obtain 
new affine structures. 

The advantage of the second type of affine structure is that it possesses MRA- 
wavelets (see |Dut2| V 

The objective of this paper is to analyze these affine structures (see section|21for 
the definition) and their " super- wavelets" , their characterization, their dimension 
function, multiplicity function and spectral function, the relation between super- 
wavelets and super-scaling functions. We follow the ideas introduced in |Dutlj . and 
the main tool will be the local trace function. 

The subspaces of (M) © ... © (R) which are invariant under all integer trans- 
lations T'^ are called shift invariant. The local trace function is a complex valued 
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map defined on M which is associated to a shift invariant space and an operator on 
P (Z) © ... (Z) f definition 14. Varying the operator, we obtain a lot of infor- 
mation about the shift invariant subspace. Taking the operator to be the identity 
yields the dimension function f definition 16. 7|) . the spectral function is obtained by 
taking the operator to be a certain canonical one-dimensional projection (definition 



The central results are theorem 14 . 41 and l4 . 51 Their power is twofold: they give a 
formula for computing the local trace in terms of normalized tight frame generators 
and also they show that this formula is independent of the choice of the normalized 
tight frame generator. Particular instances of these theorems give fundamental 
results about wavelets and shift invariant subspaces: 

(i) Take the shift invariant subspace to be (M) ® ... © (R) and the conse- 
quence is the theorem that characterizes normalized tight frames generated 
by translations: theorem 14.91 This in turn gives a characterization result 
for super- wavelets, theorem 15.41 

(ii) Take the shift invariant subspace to be Vb, the core space of the GMRA 



of a wavelet and take the operator to be /, then the result is the equality 
between the dimension function and the multiplicity function (proposition 



(iii) For Vb, take the operator to be a canonical one-dimensional projection and 
the result is the Gripenberg- Weiss formula of corollarv l6.6l 

Also, using the spectral function we are able to give a lower bound estimate on the 
dimension function (corollary I6.13|l which shows not only that some of our affine 
structures do not have MRA wavelets, but also we can give a lower bound on the 
number of scaling functions any such wavelet needs. 

Another interesting application, which follows from the characterization theorem 
15.41 is a refinement of an oversampling result of Chui and Shi |CSj . which states 
that when ^ is a NTF wavelet for N = 2, then rj := l/pij;{-/p) is also a NTF 
wavelet, p being an odd number. We prove in theorem 15 .81 that much more is true: 
f] is part of a super-wavelet fj := {i], rf) in {Wf and there is also a converse: if 
rf is a NTF super-wavelet then is a wavelet. Moreover, going from the wavelet V' 
to the super-wavelet fj preserves the orthogonality. 

Some notations: 

T is the unit circle in C. We will often identify it with the interval [— tt, tt) and 
the functions on T with functions on [— tt, tt) and with 27r-periodic functions on R. 
The identification is done via z = e~*^. 

The Fourier transform is given by 



If T in an operator on LP' (M) we denote by T its conjugate by the Fourier transform 
f / = Tf, for / e L2 (R). On (r) ^ ... ^ 2,2 (j^^ ^-^e Fourier transform is done 
componentwise. 

For / e (T), denote by 7ro(/) the operator on (k) defined by 



El- 



ion . 




Mf)^ = fV, if e (T)); 



note that Tq — t:o{z) (here z indicates the identity function on T). 
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2. Permutative wavelet representations 

We define here the affine structure on (M) ... © (R) that we wiU work 
with. This structure contains the ones used in HL and |Dut2| which we mentioned 
in the introduction. 

Let (7 be a permutation of the set {1, n} and Zfe = e^*^'-' G T, G [— tt, tt), with 
the property that — Zcr(k) ioi- all k G {1, ...,n}. We denote by Z := (zi, ...,z„). 
We define a wavelet representation on the Hilbert space 

(M)" =L^ (R)®...®L2 (K). 
" . " 

n times 

For / G (T) define the operator 7r^.z(/) on by: 

7r<T,z(/)(¥'i, = {T^o{f{zzi))ipi,...,T:o{f{zZn))ipn), 

for ipi G (M), (i G {1, which means that the Fourier transform of this 

operator is 

7fCT,z(/)(<^l, •■•,</?«) = (/(■ +Oi)ipi,...,f{- +en)(Pn)- 

Define the dilation Ua-,z on (M)" by 

Ua.zi'i^l, ■■■,Vn) = (Uoifa^i), ■■■,UQ^p^(^n))■ 
^hen (L^ (M)" , TTo-.z, ?7o-.z) is a normal wavelet representation, i.e. it satisfies the 
following conditions: 

(i) Ua^z is unitary; 

(ii) TTa^z is a unital representation of the C*-algebra L°° (T); 

(iii) U^,zna,z{f)U-}, = TT„,z{f{z^)) for aU / G (T); 

(iv) For every uniformly bounded sequence {fn)neN in L°° (T) which converges 
a.e to a function / G (T), one has that TTa.zifn) converges to TTa.zif) 
in the strong operator topology. 

Checking these properties requires just some simple computations. The translation 
(or shift) of this wavelet representation is 

Ta,z ■= T^a.ziz), 

where z indicates the identity function on T, z i— > z. Also note that z ~ '^<y.z{z^). 

Remark 2.1. Each permutation a can be decomposed into a finite product of 
disjoint cycles. Since zf^ = z„{^i) for all this implies that all Zi have a finite orbit 
under the map z i— > z^, i.e. they are cycles, as defined in the introduction. 

It may be that the length of the cycle given by Zi is shorter then the length of 
the corresponding cycle of a (for example, when all Zi — 1 and cr is a product of 
disjoint transpositions, then the length of the cycle of Zj is 1, but the cycles of a 
have length 2). However when these lengths coincide we obtain the representations 
defined in |Dut2| . 

When Zi = 1 for all i and a is the identity we have the amplification of the 
standard representation on (M), amplification which was studied in HL . 

For the rest of the paper we will consider a fixed permutation a and the points 
zi, Z2, z„, and we will omit the subscripts. 
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3. Shift invariant subspaces 

In this section we present some structural theorem for subspaces oi (M)" which 
are invariant under the integer translations T*^. 

Definition 3.1. A subspace V of (jr)" called shift invariant (SI) ii T'^V C V 
for aU k&Z. 

If ^ is a subset of L'^ (K)", we denote by S{A) the shift invariant subspace 
generated by A. 

Definition 3.2. Let F be a shift invariant subspace of (R"). A subset $ of ^ 
is called a normalized tight frame generator (or NTF generator) for V if 

{Tk^ I A: e Z, e $} 

is a NTF for V. 

We use also the notation S{(p) S{{ip}). ip is called a quasi-orthogonal gener- 
ator for S{ip) if 

{Tk'p \keZ} 
is a NTF for Sitp) and for all ^ e K, 

Per|^|2(0 :=^|^|2(e-f-2fc7r)e{0,l}. 

fcGZ 

(actually, the second condition is a consequence of the first but we include it any- 
way). 

We will need the "fiberization" techniques based on the range function of Helson 
[H] , which were used in the classical case in the work of de Boor, DeVore and 
Ron [EHRll and further developed by Ron and Shen [K^ - [RE4| . Bownik jBol| and 
others. 

For If = {cfi, ipn) G L'^ (I^)", the map Tperffi assigns to each point ^ e K a 
vector in P (Z)", called the fiber oi (p at 

Tper{(pi, fn){0 = {(M^ ~ + 2/c^))feeZ, {Vn{^ " + 2fc7r))fcez) , (C ^ , 

for aU (</7i,...,(^„) W"- 

This map transforms the translations into multiplications by scalars in each fiber: 

Proposition 3.3. For all f e L°° (T) and {(pi, ipn) e R, 

rpe.(^(/)(^l,...,V'„))(0 =/(O^P-(^l,...,^n)(0, (eeK). 

In particular, for all fc € Z, 

Note also the periodicity property of Tper (which justifies the subscript "per"): 

TpeM^ + 2sn) = A(s)*(rpe./(e)), (e e K, s e Z), 

where, for s G Z, A(s) is the shift on P (Z)": 

{\{s)a){k,i) = a{k - s,i), {keZ,ie {1, ■■■,N}). 

For a shift invariant subspace the fibers TperfiS.) at a fixed point ^ form a 
subspace of P (Z)", which will be denoted by Jper(0- Thus Jper will be a bundle 
map which assigns to each point ^ S M a subspace of P (Z)". This is the idea 
behind the range function. 
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Definition 3.4. A range function is a measurable mapping 

J : [— 7r,7r) — > { closed subspaces of (Z)"}- 

Measurable means weakly operator measurable, i.e., ^ i— > (^Pj^^-^a \ 6) is measurable 
for any choice of vectors a,b ^ P (Z)". 

A periodic range function is a measurable function 

Jper : M ^ { closed subspaces of P (Z)"}, 

with the periodicity property: 

Jper(^ + 2kn) = X{ky ( Jper(O) , [ke'L.^e R). 

Sometimes we will use the same letter to denote the subspace JperiO and the 
projection onto Jper{C)- In terms of projections, the periodicity can be written as: 

Jper(C + 2fc7r) = A(fc)* Jper(OA(fc), (fc G Z, ^ G M). 

There is a one-to-one correspondence between shift invariant subspaces and range 
functions. This correspondence is made precise in the next theorem due to Helson 
[H) . The proof given in |Bol| . proposition 1.5, carries over here, without any 
significant modification. 

Theorem 3.5. A closed subspace V of L'^ (R") is shift invariant if and only if 
V = {feL^ (R)" I Tp^rfiO e JpeAO for a.e. ^ G R"}, 

for some measurable periodic range function Jper ■ The correspondence between V 
and Jper is bijective under the convention that range functions are identified if they 
are equal a.e. Furthermore, ifV — S{A) for some countable A C (R"), then 

JperiO = spaIi{7^er'^(C) | <y5 G ^}, for a.c. ^ e R. 

Another fundamental fact about the range function is that it transforms the 
frame property into a local property: more precisely, a family of vectors generate 
a normalized tight frame by translations if and only if the fibers form normalized 
tight frames at each point in R. This is given in the next theorem, for a proof look 
in |Bol| . theorem 2.3. 

Theorem 3.6. Let V be a SI subspace of (R)", Jper 'its periodic range function 
and $ a countable subset of V . {T^^p | fc G Z, G $} is a frame with constants A 
and B for V (Bessel family with constant B) if and only if {Tper^p{£,) \ v & ^} is 
a frame with constants A and B for JperiC) (Bessel sequence with constant B) for 
almost every ^ G R. 

Another approach, with a more representation-theoretic flavor, is proposed in 
|BMM| ■ [BM| ■ [B] . The group Z has a unitary representation on a shift invariant 
space and the abstract harmonic analysis arguments given in the papers by Baggett 
et al work here with minor changes and yield the following results. 

Theorem 3.7. // Vq is a shift invariant subspace of (R)" then there exists a 
multiplicity function m : T {0, 1, ...,oo} such that: 

(i) If 

5, :={CGT|m(e)>j}, 
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then there is a unitary 

oc 

i=i 

which intertwines the representations 

(3.1) JT'' = pkJ, (fc e Z), 

where pk is multiplication by the function ^ e^'^^^ on each component: 
Pk{h,h, ■■■) = (e-*/i, e-*/2, ...), (fc e Z). 

(ii) If 

cb^■.= J-\xs,), (ze {1,2...}), 
(where xSi stands for the element of (BL'^{Sj) whose only nonzero compo- 
nent is xSi *^ L'^{Si)), then 

(3.2) {(pilie {1,2, ...}} is a NTF generator for Vb, 



(3.3) (i^j,fc,ZeZ), 
and 

(3.4) (w.(OIW.(0) = { ^o;' ff -^l (ee[-^,T). 

(iii) Vb can be decomposed as the orthogonal sum of the shift invariant subspaces 
generated by the quasi- orthogonal generators (j)i: 

oc 

(3.5) ^0=0^(00- 

1=1 

Proof. We follow the arguments given in |BMM| which use the spectral theory of 
Stone and Mackey for the commutative group Z that has as dual Z = T. 

For the representation of Z by translations on Vq, there is a unique projection- 
valued measure p on T for which 

T%, = / e-^^«dp(0, (fceZ). 

The projection-valued measure p is completely determined by a measure class 
[/i] on T and a measurable multiplicity function m : T ^ {0, 1, oo}. 

As in BMM proposition 1.2, we will show that the measure class [p] is absolutely 
continuous with respect to the Haar measure on T. 

Indeed, the Fourier transform establishes an equivalence between the translation 
T on (ffi.)" and its Fourier version 

f '=(/!, /„) = (e-^'=(-+«')/,).6{i,...,„}, {k e Z). 
Define W : (K)" ^ (R)" by 

Wifu...Jn)^if^{■-0^))^e{l,...,n}■ 

This intertwines the given representation with the one defined by 

T'^ifl, fn) = (e~*''7i)i6{l,...,n}- 
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But M can be decomposed as 

M = Ufcez[(2fc-l)^,(2fc + l)7r), 

thus, the representation of Z on the entire space (R)" is equivalent to n times 
an infinite multiple of the regular representation of Z (which is multiplication by 
e-'^^ on L^{T). 

Therefore, the projection- valued measure associated to the representation of Z on 
the entire space (M)" is equivalent to the Haar measure on T, so the projection- 
valued measure associated to any subrepresentation of this representation must 
have a measure which is absolutely continuous w.r.t. the Haar measure. 

For (|3.1|l . (|3.2() and H3.3() the argument is the one used on jBMj lemma 1.1. 

1)3.4(1 can be proved as in lemma 1.2 from |BMj . It follows from the fact that J 
is unitary and intertwining so, for k € Z^i, j G {1, 2, ...}, 

(TV. 10,) = / e-^'^^xsAjd^, 
Jt 

which implies 

Jt Jt 
H3.5|l is just a rephrasing of the previous statements. □ 



4. The local trace function 

We are able now to define the local trace function. The proofs from |Dutl| work 
here in most of the instances and, when a result is presented without a proof, the 
reader can look in |Uutl| . 

Each shift invariant subspace V, can be seen as a bundle of vector subspaces of 
l"^ (Z)" by means of a range function Jper- The local trace function associated to 
V and an operator T on P (Z)" is the function which associates to each point ^ 
the value of the canonical trace of the restriction of T to the subspace Jper(C)- We 
recall that the trace of a positive or trace-class operator T on some Hilbert space 
H, can be computed by 

Trace(T) = ^ (T/, | /,) , 

where {fi)i^i is any normalized tight frame of the Hilbert space H (see jPutlp . 

Definition 4.1. Let F be a SI subspace of (K"), T a positive (or trace-class) 
operator on P (Z)" and let Jper be the range function associated to V. We define 
the local trace function associated to V and T as the map from M to [0, oo] (or C) 
given by the formula 

Tv,TiO = Trace {TJp„{0) , e K). 

We define the restricted local trace function associated to V and a vector / in 
I' (Z)" by 

rvjiO = Trace (PfJpeAO) i= rv,P, (0), (C e M"), 
where Pf is the operator on P (Z)" defined by Pf{v) = {v \ f) /. 
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Remark 4.2. Even though, TJper{0 is not necessarily positive when T is, when 
computing the trace, we can consider instead oi TJper{^) the operator 
JperiCjT JperiC) ^ O'" <^^n usc nomiaHzed tight frames just for the subspace JperiS.) 
(not for the entire P (Z)"), in any case the expression Trace(rjper(C)) makes sense 
and is a positive number. 

When T is trace-class, since the trace class operators form an ideal, TJperiO is 
also trace class so Trace(T Jper(C)) 

Here are some ways to compute the local trace function. 

Proposition 4.3. For all f eP (Z"), 

TvjiO = WJperimW, ( for a.e. ^ G M). 

The next result is central for the paper. As we have mentioned in the intro- 
duction, this theorem can be used for two purposes: it gives a method to compute 
the local trace function and secondly, and maybe even more important, it shows 
that the formula on the right does not depend on the choice of the normalized 
tight frame. So one can look at the local trace function as an invariant for the SI 
subspace, and the theorem is an index theorem. 

Theorem 4.4. Let V be a SI subspace of (M)" and ^ C V a NTF generator 
for V. Then for every positive (or trace-class) operator T on (Z)" and any 

(4.1) tv,t{0 = E {TTp,MO I Tp.MO) , ( for a.e. ^ e R); 

(4.2) TvjiO = E I I '^peMO) I', for a.e. C e M). 

ye* 

The converse is also true: if a family of vectors can be used to compute the local 
trace then they form a NTF generator for the SI space. Moreover it is sufficient that 
this is satisfied only for some very particular vectors (see (iii) in the next result): 

Theorem 4.5. Let V be a SI subspace of L^ (K)", Jper its periodic range function 
and $ a countable subset of L^ (M"). Then following affirmations are equivalent: 

(i) (j) G V and (j) is a NTF generator for V ; 
(n) For every f eP (Z)" 

(4.3) E I (/ I TpeMO) I' = II^P-(0(/)ll', for a.e. C e R 

(iii) For every 7^ fc G Z, /, j G {1, n} and a G {0, 1, i = 

(4.4) +a'?j(e + 2Ml' = \\JperiOiSoi+aSkj)f, for a.e. C e M. 

The local trace function is related in a very simple way to another important 
invariant, the dual Gramian. This was introduced by A. Ron and Z. Shen and suc- 
cessfully used for the analysis of the structure of SI spaces and of frames generated 
by translations in |RSlj - |R,S4] . 
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Definition 4.6. Let y be a shift invariant subspacc of (M)" and <&, a NTF 
generator for V . The dual Gramian of $ is the function G* which assigns to each 
point ^ e R the matrix 

GtjjiO = E ^^(^ - + 2fc'r)^,(e - + 21^), {k, leZ,i,je {1, n}). 

We can recuperate the Gramian from the local trace function just by computing 
it at some rank-one operators. 

Proposition 4.7. If V is a SI space, Jp^r its range Junction, and $ a NTF gen- 
erator for it then, for all k,l £ G {l,...,n}, 

where Pki,ij is the rank-one operator defined by 

Proof. The first part is a consequence of theorem 14.41 For the last equality we use 
the following argument: if {fq)qeQ is an orthonormal basis for Jper(C) then 

TTace{Pki,ijJper{0) = X! (-Pfe^'j/? I /?) 

= ^{fq\ JperiO^kt) {JperiO^lj I fq) 
= {JperiO^lj I JperiO^ki) ■ 

□ 

Having these it is clear that the dual Gramian is an invariant and in fact it is a 
complete one (see |Dutl| . theorem 4.6) for the proof: 

Theorem 4.8. Let V he a SI subspacc of (M)", $i a NTF generator for V 
and $2 countable family of vectors from (K)". The following affirmations are 
equivalent: 

(i) <i>2 C and $2 is a NTF generator for V ; 

(ii) The dual Gramians G*^ and G*^ are equal almost everywhere. 

If we particularize this result to the case when the shift invariant subspace is 
(M)", then since the range function for it is constant P (Z)" at each point, 

we obtain the following characterization of families that generate normalized tight 

frames by translations: 

Theorem 4.9. Let ^ be a countable subset of L"^ (R)". The following affirmations 
are equivalent: 

(i) <^> is a NTF generator for L^ (K)"; 

(ii) 

E - - + 2fc7r) = S.,jSk, {^eR,ie {1, ...,n},k e Z). 

Next, we present some elementary properties of the local trace function. The 
proofs of these results are simple and can be found in [Dutlj for the classical case. 
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Proposition 4.10. [Periodicity] Let V be a SI subspace, T a positive or trace- 
class operator on P (Z)", f ^ (Z)". Then, for fc S Z 

(4.5) ry,T(C + 2fc7r) = Tv.x(k)T\{k)' (0: ( for a.e. ^ G R); 



(4.6) TK/(C + 2fc7r) =Ty,A(fc)/(e), {for a.e. G R), 

where A(fc) is the shift on P (Z)", i.e. the unitary operator on P (Z)" defined by 
(A(fc)a)(/, i) = a{l - fc, i) for allleZ,ie {1, n}, a € (^^"^ 

Proposition 4.11. [Additivity] Suppose (Fi)ig/ are mutually orthogonal SI sub- 
spaces (I countable) and let V — ®ieiVi. Then, for every positive or trace-class 
operator T on P (Z)" and every f (z l^ (Z)"; 



(4.7) Ty^T — ''^^TVi^T, a-e. on R; 

iei 

(4.8) '''V,/ = TVi,/, a.e. on R. 

Proposition 4.12. [Monotony and injectivity] Let V, W be SI subspaces. 

(i) V d W iff TV < Tw,T 0,-e. for all positive operators T iff Tyj < twj a.e. 
for all f eP (Z)". 

(ii) V = W iff TV J = Tw,f a.e for all f (Z)" . 

Theorem 4.13. Let (V})jgN be an increasing sequence of SI subspaces of L^ (R")? 



^ U 

T a positive (trace-class) operator on P (Z)" and f E P (Z)". Then, for a.e. ^ G R, 
"TVj.TiO increases (converges) to tv,t{0 '^^id TVjj{£,) increases to Tv.fiX). 

We will use the local trace function for the analysis of wavelets. For this reason 
we have to see how it behaves under dilation. 

Theorem 4.14. Let V be a shift invariant subspace of L^ (R)" and T a positive 
or trace-class operator on P (Z)". For I G {0, A'^ — 1} define the operators Di on 
P (Z)" by 

Aw(fc z) - I ^"^ '2kn = 2lTr -\- Oi- Ne„-^i) -^2N SIT, with s e Z, 

[ 0, otherwise, 

and the unitary operator S on P (Z)", by 

S{vi, ...,Vn) = (V<T-1(1), ...,V„-l(n)), {Vi, ...,Vn) E f (Z)"). 

Then U~^V is shift invariant and its local trace function is given by 

Tu-'v,t[0 = 2^ Tv,s'D^TDiS I — ;t — I , (4 e R)- 
1=0 ^ ^ 

Proof. By theorem 13.71 V can be decomposed into an orthogonal sum of singly- 
generated SI spaces. Also the local trace is additive and U is unitary so we may 
assume that V — S{ip) with ip a quasi-orthogonal generator. 
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Due to the commutation relation UTU ^ ~ , the vectors 

{ip^i^r :={/-ir>,|re{o,...,iV-i}} 

form a NTF generator for U^^V. Their Fourier transforms are 

VVW / ie{l,...,n} 

Then a short computation shows that, if we consider the vectors / G {0, A^- 
1} in (Z)", 

fa-^i)[- N )' " N e — - C^^-l(,) + 27rZ, 

1 0, otherwise, 



then 

Tper^-iAO = ^ ^ e— ('^ e {0, - 1}). 



N-1 

The matrix 

= e " 

/ ;,re{o,...,Af-i} 

is unitary, so the vectors 1 1 G {0, A^— 1}} span the same subspace of P (Z)" 

as the vectors {Iperf-i^riO | r e {0, iV — 1}}, namely Ju-^viO- 
Note that 

D/ are isometrics and S is unitary. In addition ip is quasi-orthogonal so ||7^erV(0 II ^ 
{0, 1} for almost every ^ G M, hence ||wi(C)ll ^ {0; 1} for all I. It is also clear that 
the vectors vi{£,) are mutually orthogonal for a fixed ^. Having all these, it follows 
that {vi{£,) 1 1 £ {0, N — 1}} is a NTF for Jjj-'-viOj therefore we can compute 
the trace with them. 

N-1 
(=0 

= 2^ {TDiSTp^rf y j;^ j \DiSTper^ ^ — 

N-1 



\S Di TDiSTperV I 7T I I ^perV 



1=0 ^ \ / \ 

which proves the formula. □ 

5. A CHARACTERIZATION OF SUPER- WAVELETS 

We apply now the local trace function to super- wavelets. Our first goal is to 
obtain a characterization of NTF super- wavelets. This is done in theorem l5.4l 

Definition 5.1. A finite subset * = {■(/'\ ■■■,ip^} of (k)" is a NTF wavelet if 
the affine system 

is a NTF for (R)". 
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The quasi-affine system is 
with the convention 

~ _ f U-^T^il), if j>0,k£Z 

'^^•'^ \ N^I'^T'^U^^i) if i < 0, fc e 1. 

There is an equivalence between affine and quasi-affine frames. This was proved 
in full generality for the (R)-case in (CSS) . Again the result will work here with 
some slight modifications in the proof. 

Theorem 5.2. Let * = t/'^} he a finite subset of L'^ (R)'' . 

(i) X(^!) is a Bessel family if and only if X'^(^!) is a Bessel family. Further- 
more, their exact frame bounds are equal. 

(ii) X('^) is an affine frame with constants A and B if and only if X'^(^) is 
a quasi-affine frame with constants A and B. 

Proof. The proof of theorem 2 in |("SSj works here word for word once we have the 
following lemma, which is the analogue of lemma 4 in 



Lemma 5.3. Let = V'l} C L'^ (M)" and let f <E L"^ (M)" be a vector such 

that /i, have all compact supports. Then 

j<o i/je* feez 



J-.^-" E E EEi(^vi^.^^)i^-o- 

Lemma 4 from |CSS| shows that this is true for the classical wavelet represen- 
tation on (K) (i.e. n ^ 1, a — identity, =0). We will show that the general 
case follows from this particular one. 

For Af > 0, j < 0, V' G ^ S Z we evaluate 



\{u-''f\i',,k)\' = \J2 / U-''f,-M^,^{x)N^/^z^T,^U~^^p,-.^,^{x)dx\' 

n „ 

<lE / \Uo''L-Mi.){x)\N^^'\T,'U,-'ij„-,^,){x)\dx\'^ 

n „ 

<lE / U^''Ff{x)N^/'T,'U^'P^{x)dx\\ 



i=l 

where 

F/ := max{|/,| |i € {1, ...,n}},P^ := maxjlV-il |« S {1, ...,n}}, 
and the last term in the previous inequality is actually 

Apply now lemma 4 of |('SSj for f ^ Ff and vJ/ ^ {P^ | V' € and with the 
previous inequality we obtain the first limit. 

The second one can be obtain by a similar argument. □ 
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Theorem 15. 21 is one of the crucial steps in obtaining a characterization for wave- 
lets. The only thing that remains to do is to use theorem 14.91 for the quasi-afhne 
system and carefully compute the Gramians to obtain the next result which gives 
the characterizing equation for NTF wavelets in (K)". 

Theorem 5.4. Let \1/ = {^i, ■■■,ipp} be a finite subset of (R)'\ The following 
affirmations are equivalent: 

(i) The affine system is a NTF for (R)"; 

(ii) The following equation hold for almost every ^ g M, hj& {1, 
For all indices i,j with Zi — Zj: 

oo 

(5-1) E E ^^."«(^'"0^<.'"0)(^"0 = %. 

ij}^^ m— — oo 

For all z, j and for all s E A^Z; 

i/iG* m>0 

Proof. According to theorem 15.21 the affine system X{'i>) is a normalized tight 
frame if and only if the quasi-affine is. This can be reformulated as 

Wm,o I m < 0, 1^ e U {^P^^r I m > 0, r e {0, iV" - 1}, e ^} 

is a NTF generator for (R)". We will use theorem 14.91 for this family. 
We have to compute the Fourier transforms: 

V / ie|l,...,n| 



and 



e{i,...,n} 



when TO > 0, r g {0, N"^ — 1}. Change to into — m and the equations of theorem 
EHIare 

E E V^--W(^"(^ - - 0, + 2fc7r)) + 

■06* m>0 

^ g-ir-(S„„(,)+7V"'(«-eO-(e<,-o)+Af"(«-e,+2fe7r))). 
r=0 

(5.3) •V^.™w(A^™(e - 0.))^,..(,)(A^"(e - + 2fc7r))) = 5^,Sl, 

i, j G {1, n}, /c e Z for almost all ^ e K. 

The rest of the proof will be just the attempt to rewrite the equation H5.3|l in 
the nicer form given in the statement of the theorem. 

We evaluate the inner sum, for fixed £,,i,j, k. Denote by 

aim) := - - 7V"(0, - 9, + 2k^), (to < 0). 
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S-^ AT™ ^ g-ir(e„™(i)+Ar'"(^-ei)-(e„™(,)+jv'"(5-e,+2fc7r))) 



r=0 



= < 1 — a(m) ' 

t 1 if e~^"(™) = 1. 

But m < so 7V^™6'o.m(i) = 6*1 mod 2tt and similarly for j so that 



Sm — 



0, if a(m) ^ mod 27r 

1, if Q:(m) = mod 2tt. 

Also note that, if for some m < 0, a(r7i) e 27rZ then a(m + 1) € 27rZ. 
Indeed, 

27V7rZ 9 Na{m) = iV6l^™(j) - NO^^^j) - N"^+'^{e, - Oj + 2A:7r) 

But 

= mod 27r 

and similarly for j therefore we see that a{m + 1) is in 27rZ also. 

Thus, if for some mo < we have a(mo) G 27rZ then all a{m) are in 27rZ for 
> m > mo, and so 5m = 1 for m > toq. 

Define 

Too := min{m < | q;(to) G 27rZ} G {— oo, —1, 0}. 

Rewrite (|5l^ : 

'06^>I' m>mo 

Of course, toq depends on i, j, k and we have to make this dependence more precise. 

If Zi = Zj and fc = then 6am.[i) = d<T"'{j) for all m G Z so a (to) = for all 
TO < and therefore mo = — oo. In this case the equation (|5.4|l is 



m— — OQ 



which after the change of variable £, — Oi £^ becomes (|5.1|l . 

If Zi Zj or k =/= then mo is finite. When Zi ^ Zj we have 0a"^[i) ^ ^^"^(j) for 
all TO, by cyclicity 9a-rn(^i^ — Ocrm(^j-) will assume only a finite number of nonzero values 
in (— 27r, 2tt) as to varies. For m close to — oo the quantity —N™{di — 6j + ikir) is 
too small to cover the distance from O^m^j^ — (^cr"^{j} to any point in 27rZ so that 
a(m) ^ 27rZ. 

We change the variable in C - 6*^-0 (i) = N"^o{^ -Oi). Then 

iV™0(^ „ 0^ + 2kTT) ^a-oO) - a(TOo), 

so (|5.4f) is equivalent to 

(5.5) E V^.™+"o(.)(iV'"(C - 0„^oi^)))HN'^{^ - e,..oU) - a(TOo))) = 0, 

m>0 

We prove that (|5.5() is equivalent to: 
For all i,j and s G Z \ A^Z: 

(5-6) E E^-'"W(^"(^-^'))^-'"0)(^'"(^-^^ +^(*'j') + 2^'^)) = 0' 

i/'e* m,>0 
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where A{i,j) = iV(0,-i(,) - ^a-iy)) " - ^j)- 

First, assume (I5.6|) holds. Fix i,j S {1, n} , k E Z with Zi ^ Zj or 7^ 0. We 
prove that -a(mo) is of the form A(cr"°(i), cr™°(j)) ^ 2s7r with s e Z \ 

We have 

a{mo) = 0,-^0 (,) - 0,^0 u) - N"">{0^ - 9, + 2kTr). 
By the definition of niQ, a{mo — 1) ^ 27rZ. But observe that 

Naima - 1) = iV(0,™„-i(.) - ^.^o-io)) - N"''>{9, - 9, + 2kTi) 
^a{m^)+A{a^»(i),a'^-{j)). 

Also A(cr™«(i),cr"«(j)) e 27rZ and a(mo) G 27rZ so Na{mo - 1) G 27rZ which 
means that a(mo — 1) G 7^^- 

We also know that Q!(mo — 1) ^ 27rZ, hence a(mo — 1) = 27r-^ with s G Z \ NIj. 
Then 

-a(mo) = -2s7r + A(a'""(z),a"°(j)). 
Now using (|5.t)|) for o-™" (i), cr"'" (j') and — s, we obtain H5.5|l . 

For the converse, suppose H5.5|l holds and fix i,j and s G Z \ NZ. Let k := 
2s7r + j) and we show that the mo associated to i, j and k is 0. 

Indeed, we compute: a(— 1) = — 2-^7r ^ 27rZ. Since a(0) = — 2s7r — j4(i, j), using 
(|5.5(l . the relation l|5.6|l is obtained. 

Equation (|5.()|) is clearly equivalent to H5.2|l after a change of variable. □ 

Remark 5.5. Note that there is some redundancy in the equations of theorem 
15.41 Indeed, the relations H5.1|l and H5.2|l for i > j follow from the corresponding 
relations for i < j after a conjugation and a change of variable. 

Also, if ii, ii form a cycle for a then, once we have equation H5.1(l for i — j — ii, 
the equation 1)5.1(1 will be also true for 12, because one needs only to change 
the variable ^ <-> N^^ for some appropriate p. 

We can apply theorem 15.41 to the particular case of the amplification of the 
classic representation on (R). Our result complements the equations obtained 
for orthogonal super- wavelets in HL , theorem 5.13. 

Corollary 5.6. Let a be the identity permutation and zi = ... = z„ = 1. Then ^ 
is a NTF wavelet if and only if the following equations are satisfied for almost all 
C G R: 

00 

J2 E MN"'O^^A^'"0=kj, {i,je{l,...,n}); 

ip^^ m— — 00 
00 

E E^»(^'"^)^j(^"(^ + 2s'^))=0' G {l,...,n},5GZ\iVZ). 

In |HL| . two normalized tight frames (ei)ig/ for H and {fi)i£i for K are called 
strongly disjoint if their direct sum (e^ © fi)i£i is a normalized tight frame for 
H (B K. Two normalized tight frame wavelets 5* = {V^i, V'l}, ^' = {ip'i, V'lI 
for (R) (note that the cardinality is the same) are strongly disjoint if their affine 
systems X{tfj) and X{'i>') are strongly disjoint. 

If we look at corollarv l5.6l we see that the fact that each component {ipi \ ip G 5*}, 
(i G {1, n}) is a NTF wavelet for (R) is equivalent to the equations with i = j 
(and we reobtain in fact the well known characterization of wavelets in (R), see 
|Bo2j ■ [Caj or |HWp . Therefore the disjointness part is covered by the equations 



LOCAL TRACE FOR SUPER- WAVELETS 



17 



with i ^ j 1 hence we can use the coronary 15.61 to produce a characterization of 
disjointness: 

Corollary 5.7. Two NTF wavelets * {■(/'\ V'^}, := V'''^} for 

LP' (K.) are strongly disjoint if and only if 

L OQ 

l — l m— — oo 
L CO 

^ J2 V^'(^'"OV^(^"(C + 2s7r)) = 0, (s e Z\iVZ). 

i = l m=0 

In the remainder of this section we present another apphcation of the character- 
ization theorem. The following oversampling result is proved in |CSj for the scaling 
factor A?^ = 2: if one has a NTF wavelet (for (k))^ then also rj := ^ip is a 
NTF wavelet, where p is any odd number. 

We will refine this result here: not only that 77 is a NTF wavelet, but it is part of 
a super- wavelet (77, ry) for (K)^, and the starting wavelet is an orthogonal one 
if and only if the super-wavelet is also orthogonal. Here is the precise formulation 
of this result: 

Let p be a positive integer which is prime with N. For ?/; G (K), define 
ry(7/.)(x) = -^J (-] , {x e R), 

rJiij) = (r^i^P), ry(V')) G [Rf , € 
Let p :— e~^'^*/P, Zk — p'", k g {1, ...,p} and a the permutation of {1, defined 
by (p^)^ = /o'^*-*'' for all i e and let Dla.z be the corresponding affine 

structure on (R)^. 

Theorem 5.8. Let * := {ipi, ...,ipp} in (r). 5» is a A^TF waDefet /or (k) 
anrf onZy if 

^(*) := {7/(V') I V e *} 
is a NTF wavelet on (M.)^ for the affine structure K^.z- 

Moreover, 4' is an orthogonal wavelet for (R) if and only if f]{'^) is an or- 
thogonal wavelet for (R)^. 

Proof. We have to check that the equations of theorem 15.41 are satisfied. 
The Fourier transform of ri{'ip) is 

ri^){0 = ${px), (eeR). 

Then, we can rewrite the equations (j5.1fl and (j5.2|l for f]{ip): 

00 

■0^^ m— — 00 

E E ^^(^^>o^(^™(pe + «)) = 0, 

i/>e* ■m>0 

for all a e {p(A^6'^-i(,) - NO^^iq) + 2s7r) | i, j e {1, s e Z \ iVZ} =: A. We 

want to see what this set is and we will show that it is equal to {2g7r | g e Z \ A^Z}. 



18 



DORIN ERVIN DUTKAY 



Note that 0i takes all the values of the form 2kn/p inside [— tt, tt) with k integer. 
Then pOi takes all values of the form 2k-K inside [— p7r,p7r) and therefore p0a-i(i) — 
pO„-i(^j) will cover all values of the form 2fc7r within (— 2p7r, 2p7r). Therefore 

A = {2tt{NI+ps)\1 C + 1 - p + 2,...,-l,0, l,...,p- l},s e Z\ A^Z}. 

Since p is prime with N it is clear that all numbers in A are of the form 2'iTq 
with 9 G Z \ NZ. 

For the converse, take q not divisible by N . Since N and p are prime we can 
write q = Nli + psi with /i, si integers. Clearly si cannot be divisible by N . Also 
we can write li = pri + I with ri e Z and I g {0, — 1}. So 

q^ Nl+p{Nri+si) 

with i G {0, ...,p — 1}, A^ri + si not divisible by and therefore 27rg is in A. 

Finally, if we change the variable p^ ^, we see that the equations are equivalent 
to the fact that * is a NTF wavelet for (R). 

The orthogonal case follows by an inspection of the norms: a normalized tight 
frame is an orthonormal basis if and only if all the norms of the vectors are 1. But 
observe that 

and therefore the last equivalence is clear. □ 

6. The spectral function and the dimension function 

Definition 6.1. A finite subset = of (R)" is called a semi- 

orthogonal wavelet if the afBne system 

{t/^rVli e Z,fc e Z, V e 

is a NTF for (R") and W^ _L Wj for j, where 

Wj = Vj+i e Vj = span{C/-^T'=?A | fc G Z, V e = C/-^'(5(^')), (j G Z). 

and {Vj)j is the GMRA associated to 

If (f> is a NTF generator for Vq then $ is called a set of scaling functions for '5. 

For a semi-orthogonal wavelet, theorem 15.21 and the orthogonality conditions 
given in the definition imply the following: 

Proposition 6.2. Let be a semi- orthogonal wavelet and (Vj)j its GMRA. Then 
the vectors 

{i^m,0 = iV-"/2[;m^ | m > 0, ^ G *} 

form a NTF generator for Vq . 

This proposition will be one of the essential ingredients for our results. In Vq we 
have two NTF generators: one given by the wavelet VE" as in proposition 16.21 and 
another given by a set of scaling functions for Vq which always exists according to 
theorem 13. 71 Computing the local trace for various operators in two ways will give 
rise to some fundamental equalities about wavelets. 

The spectral function was introduced for the classical case on (R) by M. 
Bownik and Z. Rzeszotnik in |BoB,zj . It contains a lot of information about the 
shift invariant subspace and it has several nice features. In some sense, the spectral 
function corresponds to the diagonal entries of the dual Gramian. 
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Definition 6.3. Let y be a shift invariant subspace of (M)". For i g {1, 
the z-th spectral function of V is 

where Pqj is the projection onto the 0-th component of the z-th vector 

iPQ,iv)ik,j) = wo.j^ofe^y , for aU j e {1, ...,n},k e Z. 

The next proposition is just a consequence of theorem 14.41 

Proposition 6.4. Let V be a shift invariant subspace and ^ be a NTF generator 
for V . Then 

'^vM) = E - ^')' e K, z e {1, ...,n}). 

Using the NTF generator given in proposition 16.21 we obtain a formula for the 
spectral functions associated to wavelets: 

Proposition 6.5. Let be a semi- orthogonal wavelet and {Vj)j its GMRA. Then 
the spectral functions of the core space Vq can be computed by: 

Combining proposition 16 . 41 and 16 . 51 we obtain a Gripenberg- Weiss- type formula: 

Corollary 6.6. //^ is a semi- orthogonal wavelet and $ is a set of scaling function 
for it, then 

The dimension function was extensively used on (K) to decide whether a 
wavelet is a MRA wavelet or not, see |HW,. It also gives information about the 
multiplicity of the wavelets r[Bl.[BM|. [BMM] . |V\^ 'l. 

Definition 6.7. Let y be a shift invariant subspace of L^ (M)" and Jp^r its range 
function. The dimension function of V is 

dimv.(0 = dim(Jpe.(e)), (e e K)- 

Since the trace of a projection is the dimension of its range, using the definition 
of the local trace we have the following equality: 

Proposition 6.8. Let V be a shift invariant subspace of L^ (M)". Then 

dimviO = TvjiO, 
where L is the identity operator on P (Z)". 

Then, with theorem 14.41 and proposition 16.41 we have 
Proposition 6.9. //$ is a NTF generator for a shift invariant subspace V then 

n n 

dtmv{o = EEEi^^i'(^-^^+2M = EPe^('^v-,.)(e), (eeM). 

i=l fegz i=l 

where 

Per(/)(e) - E + 2*^^)' (e e K, / e (K)). 

fcez 
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We can choose the scahng functions to be as in theorem 13.71 and the conse- 
quence is that the dimension function and the multiphcity function are the same 
(a generahzation of the resuh of E. Weber | Webj ) . 

Theorem 6.10. // 5* is a semi- orthogonal wavelet and Vq the core space of the 
associated GMRA. Then the multiplicity function and the dimension function ofVo 
are equal and they are given by 

n 

i/>e* m>i i=i fcez 
for almost all f € [— tt, tt). 

Next we want to estabhsh some dilation formulas for the spectral function and 
for the dimension function. We use theorem 14. 141 

Proposition 6.11. Let V be a SI subspace of (M)". Then for all j e {1, ...,n}, 

<Tu-^vAO = <^v,a-Kn i^- — ' ^ ^^' j , (e e M), 

dimjj-iviO = E ^^"^^ ( — ]v — ) ' ^ 
1=0 ^ ^ 

Proof With theorem Ol we have to compute for I e {0, TV - 1}, S*D*iPo^jDiS 
is the rank-one operator given by the vector S*D'[eQ^j. We have for k E Z,i E 
{0,...,iV-l}: 

S*D*eo.j{k,i) = {D*eo,j){k,a{i)) = eo,j{2lTT + 0„^i-, - N9, + 2kNn), 
and this shows that 

S*DreQ = I %.--Mj) if l = h 

[ 0, otherwise, 

where Ij E {0, ...,N — 1} and kj E Z are the only such numbers that verify the 
equation 2ljTT + dj — N9„-i(^j-f + 2kjNTT = 0. Therefore 

S*D*P DS = I ^''i^'^-^ii) " ^(^■j)Po,a-^])Mko)* if I = Ij, 
\ 0, otherwise. 

Using theorem 14 ■ 1 41 and then the periodicity property given in proposition l4.10l we 
obtain: 



f^ + 2l^\ 



f^ + 2lj7r 
= crvo,a-i(i) ( ]^ 1" 

and, with the definition of kj and Ij the equation follows. 

The equation for the dimension functions follows from the fact that S* IDiS = 
/ for alU e {0,...,7V- 1}. □ 

Proposition 6.12. Let {Vj)j be a GMRA and $ a NTF generator for Vq. Then, 
for almost every ^ G R and all i E {1, n}, 

lim y ( J- 
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Proof. With the monotone convergence theorem 14.131 we have that cru-^^v,i{0 
converges pointwise a.e. to cr2^2(R')n j(^) = 1. With proDOsition 16 . 1 ll by induction 
we get: 

iV™ J 
Changing ^ — 0^ to ^ we obtain 

But (7 is a finite permutation so cr^ — id for some p. Thus, for the subsequence mp, 
a~"^P{i) = i for all m, hence 

Jim^^l^.p =1, 

and if we apply this to ^, ^/N,...,^/NP^^ we obtain the desired limit. □ 

The limit given in proposition 16.121 enables us to give a lower bound for the 
dimension function of a wavelet. We can conclude that, if there are two distinct 
indices i,j such that Zi — Zj then the dimension function of any wavelet is strictly 
bigger then 1 at some points so there are no MRA wavelets, and for any wavelet 
one needs at least 2 scaling functions. If there are three indices for which the points 
Zi are the same, then any wavelet needs at least three scaling functions, and so on. 
This result generalizes and refines proposition 5.16 in jHL| . 

Corollary 6.13. // 5* is a semi- orthogonal wavelet then, for all G [— 7r,7r), 

limsup£'>i((^) > card{i G {1, n} \ 9i — uq}. 

In particular, if there are two distinct indices i ^ j with Oi — 9j then there are no 
MRA wavelets. 

Proof. Let $ be a set of scaling functions associated to ^' and compute I?* at 
^/TV™ + ao using proposition 16. 91 Then take the limsup and use proposition 16. 121 
For a MRA wavelet is constant 1, and the second statement follows by 
contradiction. □ 
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